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Abstract 

Using 5.63 pb^^ of data accumulated at the il^{2S) resonance with the CLEO III and CLEO-c 
detectors corresponding to 3.08 million il>{2S) decays, a search is performed for the decay ip{2S) 
r/cTT"*" 7r~7r° to test a theoretical prediction based upon the assumption that the cc pair in the 
ip{2S) does not annihilate directly into three gluons but rather survives before annihilating. No 
signal is observed, and a combined upper limit from six r/c decay modes is determined to be 
B{'ip{2S) rjcTT^TT^TT^) < 1.0 X 10^^ at 90% C.L. This upper limit is about an order of magnitude 
below the theoretical expectation. 
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In perturbative QCD the charmonium states J/ip and ipi^S) are nonrelativistic bound 
states of a charm and an anticharm quark and it is predicted that the decays of these states 
are dominated by the annihilation of the charm and anticharm quark into three gluons. 
The partial width for the decays into an exclusive hadronic state h is then expected to be 
proportional to the square of the cc wave function overlap at zero quark separation, which 
is well determined from the leptonic width [1]. 

Since the strong coupling constant, as, is not very different at the J/ip and ip{2S) masses, 
it is expected that for any state h the J/ip and ip{2S) branching fractions are related by ^ 



where B denotes a branching fraction, and the leptonic branching fractions are taken from the 
Particle Data Group (PDG) [l|. This relation is sometimes called "the 12% rule". Modest 
deviations from the rule are expected Although the rule works well for some specific 
decay modes, isospin conserving ip{2S) decays to two-body final states consisting of one 
vector and one pseudsoscalar meson exhibit strong suppression: pvr is suppressed by a factor 
of seventy compared to the expectations of the rule (the so-called pvr puzzle) [H, 0, 0, 0] . Also, 
vector-tensor channels such as pa2(1320), and ii'*(892)K'|(1430) are significantly suppressed 
[l|, 0]. Another issue is the hadronic excess in ip{2S) decays: the inclusive hadronic decay 
rate of ipi'^S) is larger than that expected from an extrapolation of the J/ijj hadronic decay 
branching fraction by 60 — 70%. A recent review [3| concludes that current theoretical 
explanations of ip{2S) decays are unsatisfactory and that more experimental measurements 
are desirable. 

A recent paper ^ suggests that "survival before annihilation" could be an important 
mechanism in ip{2S) decays. This model proposes that the cc pair in the ip{2S) does not 
annihilate directly into three gluons but rather "survives" before annihilating, i.e., the cc 
pair decays by first emitting two or three nonperturbative gluons before annihilating into 
three or two perturbative gluons. This model, it is claimed, can solve the problem of the 
apparent hadronic excess in ip{2S) decays as well as the pvr puzzle. One important prediction 
of the model is that the ip{2S) — >■ rjcTT^n'TT^ channel would be a significant decay with a 
branching fraction of 1% or larger. 

We search for the decay ^'(25') r/cTr+Tr^vr^ using the six decay modes of the r/c listed in 
Table [H The t]c decay modes selected amount to about 10.6% of the total r]c decay rate [if. 

The data sample used in this analysis was obtained at the ip{2S) in e"'"e~ collisions pro- 
duced by the Cornell Electron Storage Ring (CESR) and acquired with the CLEO detector. 
The CLEO III detector ^ features a solid angle coverage for charged and neutral particles 
of 93%. The charged particle tracking system, operating in a 1.0 T magnetic field along the 
beam axis, achieves a momentum resolution of ~0.6% at p = 1 GeV/c. The calorimeter 
attains a photon energy resolution of 2.2% at ii^^ = 1 GeV and 5% at 100 MeV. Two particle 
identification systems, one based on energy loss {dE/dx) in the drift chamber and the other 
a ring- imaging Cherenkov (RICH) detector, are used together to separate from vr^. The 
combined dE/dx-HlG]^ particle identification procedure has efficiencies exceeding 90% and 
misidentification rates below 5% for both vr^ and for momenta below 2.5 GeV/c. 

Half of the ip{2S) data were taken after a transition to the CLEO-c \1Q\ detector con- 
figuration, in which the CLEO silicon small radius tracking detector was replaced with a 
six-layer all-stereo drift chamber. The two detector configurations also correspond to differ- 
ent accelerator lattices: the former with a single wiggler magnet and a center-of-mass energy 
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spread of 1.5 MeV, the latter (CESR-c [13]) with six wiggler magnets and an energy spread 
of 2.3 MeV. 

The integrated luminosity (£) of the datasets was measured using e~^e~ , 77, and 
final states Event counts were normalized with a Monte Carlo (MC) simulation based 
on the Babayaga [13] event generator combined with GEANT-based [3] detector modeling. 
The data consist of £=5.63 pb~^ on the peak of the ip{2S) at ^/s = 3.686 GeV (2.74 pb~^ 
for CLEO III, 2.89 pb'^ for CLEO-c). 

Standard requirements [14] are used to select charged particles reconstructed in the track- 
ing system and photon candidates in the Csl calorimeter. We require tracks of charged 
particles to have momenta p > 100 MeV/c and to satisfy | cos^^| < 0.90, where 6 is the 
polar angle with respect to the e+ direction. Each photon candidate satisfies > 30 MeV 
and is more than 8 cm away from the projections of tracks into the calorimeter. Particle 
identification from dE/dx and the RICH detector is used for all charged particle candidates. 
Pions and kaons must be positively and uniquely identified, i.e., pion candidates must not 
satisfy kaon or proton selection criteria, and kaon candidates obey similar requirements. 

The invariant mass of the decay products from the following particles must lie within lim- 
its determined from MC studies: 7r° (120 < M^^ < 150 MeV), r/ (500 < M^^ < 580 MeV), 
7] (530 < M^+^-^o < 565 MeV). 

For 7r° — s> 77 and rj ^ 'j'y candidates in events with more than two photons, the com- 
bination giving a mass closest to the known vr*^ or r] mass is chosen, and a kinematically 
constrained fit to the known parent mass is made. Fake vr*^ and t] mesons are suppressed 
by requiring that each electromagnetic shower profile be consistent with that of a photon. 
For 1] 7r+7r^7r° the vr^ is selected as described above, and then combined with all possible 
combinations of two oppositely charged pions choosing the combination that is closest to 
the rj mass. A kinematically constrained fit is not used for this mode. 

Events having final state particles consistent with one of the six r]c decay channels and 
additionally a tt+tt^tt'^ combination are selected. Energy and momentum conservation re- 
quirements are then imposed on the event using the summed vector momentum J2 Pi ^"^^ 
scalar energy ii^vis- These requirements are based on the experimental resolutions as deter- 
mined by Monte Carlo for each of the final states. The scaled momentum | I]Pi|/-^c.m. is 
required to be consistent with zero and the scaled energy E^i^/E era. is required to be con- 
sistent with unity. The experimental resolutions are less than 1% in scaled energy and 2% 
in scaled momentum. In order to veto the final states ip{2S) —>■ XJ/ip {X = 7r"'"7r~, tt^tt", 
or Tj) events are rejected in which the mass of any of the following falls within the range 
3.05 < m < 3.15 GeV: the two highest momentum oppositely charged tracks, the recoil mass 
against the two lowest momentum oppositely charged tracks, or the mass recoiling against 
the TT^TT^ or rj. 

The efficiency, e, for each final state is the weighted average obtained from MC simu- 
lations for both detector configurations. No initial state radiation is included in the 
Monte Carlo, but final state radiation is accounted for. The efficiencies in Table |T] include 
the branching fractions of the 77 [l|. 

After the selection of events consistent with the six exclusive ip{2S) decay modes, we 
search for r]c production in these events by examining the invariant mass of combinations of 
particles consistent with an r]c decay. In all events there are multiple r]c combinations and 
to be conservative we have taken only one combination per event choosing the one nearest 
the Tjc mass. Figure [1] is the scaled energy distribution for each of the six ip{2S) exclusive 
decays after imposing momentum conservation (| I]pi|/-Ec.m. < 0.02) and for events with 



4 



an rjc candidate mass greater than 2.7 GeV, which includes the nominal rjc mass region. 
There is clear evidence of exclusive ip{2S) production of the final states under study from 
the accumulation of events near unity. After requiring the scaled energy to be in the range 
(0.98 — 1.02), the invariant mass distributions shown in Fig. [2] are analyzed for rjc produc- 
tion. A combined mass distribution for all six modes is shown in Fig. [31 The histograms 
show the expected signal from Monte Carlo normalized to the branching ratio prediction 
of B{ip{2S) — >■ rjc'K^'K'TT^ = 1% . There is no evidence for rjc production in any of the six 
individual decay modes or in the combined distribution. 

Table [T] shows the number of events in the r]c mass region (2.91 — 3.05 GeV), which 
corresponds to ±3 standard deviations of the expected rjc mass resolution, the efficiency and 
the rjc branching ratio for each of the decay channels. No evidence is found for ip{2S) 
rjcT^' 
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FIG. 1: The scaled total energy distribution for each ip{2S) mode for events with the candidate rj, 
decay having an invariant mass greater than 2.7 GeV. a) rjc — > K~^K~-k^ . b) rjc ■ 
c) Vc ~^ rjTT^TT" ,rj — > vr^vr^vr'^. d) K^K^tt^tt^ . e) 77c vr+vr^Tr+Tr^. f) rj, 



rjTT^TT , rj ^ 
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TABLE I: No evidence is found for ^{2S) — > r/cVr+Tr vr''. Number of events in the rjc mass region, 
detection efficiency, and rjc branching ratio 



decay mode 


T^events 


efficiency (%) 


SiVc) (%) 


rjc K+K'TT^ 


1 


3.08 


1.2 ±0.2 


rjc rjTT+Tr~,rj 77 





2.76 


1.3 ±0.5 


rjc rj7r~^TT~,rj — > tt~^it~tt^ 





0.80 


0.7 ±0.3 


rjc K~^K~Tr~^TT~ 


7 


3.07 


1.5 ±0.6 


rjc vr+vr^vr+TT^ 


6 


4.06 


1.2 ±0.3 


rjc K-TT+K^ 


4 


1.55 


4.7 ±0.8 


all modes combined 


18 


2.31 


10.6 ± 1.2 
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FIG. 2: The invariant mass distribution for each r]c mode. Histogram: Monte Carlo (y axis on left), 
normalized to B{il){2S) rjcir^ ir^ tt^) = 1%; filled circles with error bar: data {y axis on right), 
a) Tjc — > K~^K~Tr^. b) r]c r/vr+vr^, rj — > 77. c) tJc — > t^tt+vt^, t] tt^tt^tt^. d) rjc K^tt^tt^ . 

e) r]c — > TT+vr^vr+vr^. f) r/c K^tt^K^. 
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FIG. 3: The invariant mass distribution for all six rjc modes combined. Histogram: Monte Carlo, 
normalized to B{ip{2S) — > T/cTr+vr^vr'') = 1%; filled circles with error bar: data. 



The following sources of systematic uncertainties have been considered. The number of 
ipi^S) decays (3.0%), trigger efficiency (1.0%), and the uncertainty associated with the res- 
olution defining a signal region in the scaled energy and resonant mass (2.0%) contribute 
identical systematic uncertainties to each channel. Other sources vary by channel, for exam- 
ple, Monte Carlo statistics (2.7% — 7.4%). The systematic uncertainty associated with the 
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charged track finding is 0.7% per track. This uncertainty adds coherently for each charged 
track in the event. The particle identification uncertainty is 0.3% per charged pion, and 
1.3% per charged kaon. This uncertainty adds coherently for each charged track in the 
event since PID is used for every track. The systematic uncertainty associated with 7r° and 
7] ^ finding is 4.4%. 

We determine a ip{2S) —* rjcir^Ti^Ti^ branching fraction upper limit at 90% CL. combining 
all ?7c decay modes using 

Nq 

B = — (2) 

(E^.-i3(r/,-^/iO)-^^(2S)' ^ ^ 

where Ns is the upper limit on the number of signal events, Si and B{ric hi) are the 
efficiency and the branching fraction for rjc decay mode i; and N^(2S) is the total number of 
ip{2S) decays. The number of 'tpi^S) decays N^(2S) was determined to be 3.08 x 10^ by the 
method described in 15|. The combined efficiency, which includes the branching fractions 
for Tjc decays, is defined as 

^ S gj • B{r]c -» hi) 
'~ EBivc^K) ' 

and is also shown in Table [H To determine the upper limit on the number of signal events, 
Ns, the data distribution in Fig. [3] is fit using a polynomial and a Gaussian with fixed signal 
shape obtained from MC. By fixing the signal amplitude and maximizing the fit likelihood 
with respect to the parameters for background, we obtain the fit likelihood as a function of 
the signal amplitude and determine the upper limit on the number of signal events to be 
< 6.7 at 90% CL. After combining the systematic uncertainties and the uncertainties in the 
rjc branching ratios in quadrature and increasing the upper limit on the number of observed 
events by one standard deviation of the combined uncertainty, the branching fraction upper 
limit is computed to be 1.0 x 10^'^. [Alternative methods of computing a 90% CL. upper 
limit using the number of events in the signal and sideband regions of the invariant mass 
give slightly lower limits (H. [l6|.] 

In summary we have searched for the decay ip{2S) tIcTt^it^tt^ using six decay modes of 
the ?7c. We have determined a combined upper limit at 90% confidence level for the branching 
fraction of the decay ipi^S) — ?7c7r+7r~7r° to be 1.0 x 10~^. This upper limit is about an 
order of magnitude below the theoretical prediction of the survival before annihilation model 
meaning that the survival of the cc pair to form an rjc is highly suppressed. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, and the Natural Sciences and Engineering Research Council of Canada. 



[1] Particle Data Group, W.-M. Yao et ai, J. Phys. G 33, 1 (2006) 

[2] W. S. Hou and A. Soni, Phys. Rev. Lett. 50, 569 (1983); G. Karl and W. Roberts, Phys. Lett. 
B 144, 263 (1984); S. J. Brodsky et al, Phys. Rev. Lett. 59, 621 (1987); M. Chaichian et 
al, Nucl. Phys. B 323, 75 (1989); S. S. Pinsky, Phys. Lett. B 236, 479 (1990); X. Q. Li et 
al, Phys. Rev. D 55, 1421 (1997); S. J. Brodsky and M. Karliner, Phys. Rev. Lett. 78, 4682 
(1997); Yu-Qi Chen and Eric Braaten, Phys. Rev. Lett. 80, 5060 (1998); M. Suzuki, Phys. 



7 



Rev. D 63, 054021 (2001); J. L. Rosner, Phys. Rev. D 64, 094002 (2001); J. P. Ma, Phys. 

Rev. D 65, 097506 (2002); M. Suzuki, Phys. Rev. D 65, 097507 (2002). 
[3] Y. F. Gu and X. H. Li, Phys. Rev. D 63, 114019 (2001). 
[4] CLEO Collaboration, N. E. Adam et al, Phys. Rev. Lett. 94, 012005 (2005). 
[5] EES Collaboration, M. Ablikim et al, Phys. Lett. B 614, 37 (2005), Phys. Lett. B 598, 149 

(2004), Phys. Rev. Lett. 93, 112002 (2004), Phys. Rev. D 72, 012002 (2005); Phys. Lett. B 

619, 247 (2005), and Phys. Rev. D 70, 012003 (2004); BES Collaboration, J. Z. Bai et al, 

Phys. Rev. D 70, 012006 (2004). 
[6] F. A. Harris, Int. J. Mod. Phys. A 20, 445 (2005). 
[7] BES Collaboration, J. Z. Bai et al, Phys. Rev. D 69, 072001 (2004). 
[8] P. Artoiscnct et al, Phys. Lett. B 628, 211 (2005). 

[9] CLEO Collaboration, Y. Kubota et al, Nucl. lustrum. Methods Phys. Res., Sect. A 320, 

66 (1992); D. Peterson et al, Nucl. lustrum. Methods Phys. Res., Sect. A 478, 142 (2002); 

M. Artuso et al, Nucl. lustrum. Methods Phys. Res., Sect. A 554, 147 (2005). 
[10] CLEO-c/CESR-c Taskforces & CLEO-c Collaboration, Cornell University LEPP Report No. 

CLNS 01/1742, 2001 (unpublished). 
[11] CLEO Collaboration, G. Crawford et al, Nucl. lustrum. Methods Phys. Res., Sect. A 345, 

429 (1994). 

[12] C. M. Carloni Calame et al, Nucl. Phys. B, Proc. Suppl. 131, 48 (2004). 
[13] Computer code GEANT 3.21, in R. Brun et al, CERN Report No. W5013, (1993) (unpub- 
lished). 

[14] CLEO Collaboration, R. A. Briere et al, Phys. Rev. Lett. 95, 062001 (2005). 
[15] CLEO Collaboration, S. B. Athar et al, Phys. Rev. D 70, 112002 (2004). 
[16] G. J. Feldman and R. D. Cousins, Phys. Rev. D 57, 3873 (1998). 



8 



